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p | ■ In coupled chaotic bistable systems such as Lorenz and Chua oscillators, 

two-phase domains corresponding to the two lobes of the strange attractor 
are formed. The dynamics of each domain is confined to one lobe and typ- 
ically exhibits one of the two types of behavior: oscillation death or nearly 
■ periodic oscillations. We elucidate the role of intrinsic broad time scales on 
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the confinement and the oscillation death. 
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Many systems in physics, chemistry, and biology can be considered as assemblies of cou- 
pled oscillators These systems can be categorized by the local dynamics of oscillators 
and the network structures between them. In reaction diffusion systems, oscillators are cou- 
pled diffusively, and the local dynamics can be bistable, excitable or oscillatory. They are 
classified according to their local dynamics and the characteristics of each class have been 
well understood ||. 

In recent years, cases where the local dynamics is chaotic have aroused a great deal of 
interest. Studies were performed using various models such as coupled map lattices J| -f§, 



and coupled oscillators HHUJ. Synchronization [|T3Hl5|] is the most important feature of two 
coupled chaotic oscillators. However N coupled chaotic oscillators with large N exhibit not 
only synchronization but also more complex behaviors such as nontrivial collective behavior 
H|7| and size instability ||. In particular, it has been revealed that from topological con- 
siderations line defects are a general property of coupled oscillators with period doubling 
bifurcations |9|], and such line defects have been observed in a recent experiment [ |T6|j . 

However there are many other chaotic oscillators which do not follow period doubling 



bifurcations [T7| and more studies on coupled systems of such oscillators are needed. Es- 



pecially we want to know whether universal properties exist and can be categorized by the 
characteristics of the local dynamics. In this paper, we investigate the dynamics of the N 
coupled chaotic bistable systems. We call systems like the Lorenz, Chua, and Duffing oscil- 



lators "chaotic bistable systems" [15j . In these systems, two bistable basins exist below some 
threshold, and chaotic orbits wandering between two regions emerge above this threshold. 

To begin with, we consider the following N diffusively coupled Lorenz oscillators with 
periodic boundary conditions, 

-J 1 = a{ui - Zi) + e(x i+1 + - 2xt) (1) 

%- = rxi - yi - XiZi + e(y i+1 + yi-i - 2y l ) 
at 

dzi 

— j = Xiyi - bzi + e(z i+1 + z^ x - 2^). 
at 



We will concentrate on the canonical parameter values a = 10, b = 8/3 and r = 28 [19 



The uncoupled Lorenz oscillator shows a chaotic trajectory consisting of two lobes. The 
trajectory circulates around a right saddle-focus a few times, moves to the left, and circulates 
around a left saddle-focus, and so on. These alternations are repeated indefinitely and 
erratically. 

One thing to note for later reference is that the Lorenz oscillator has broad time scales 
l20fl . Characteristics of a chaotic attractor can be well resolved by unstable periodic or- 



bits(UPO) |2T[ and UPOs of Lorenz oscillators show various periods pO| . Another method 



to examine the broad time scales is to calculate the instantaneous frequency uj{t). For Lorenz 
systems, the phase variable <ft can be defined by tan _1 (z — zq)/(A — Aq) where Zq — r — 1, 



A = y/x' 2 + y 2 , Aq = y2b(r — 1) and then uj{t) = 0(t) is given by 

u\ - A(A-Aq)z- (z - z )(xx + yy) 
U ® ~ M(z - z ) 2 + (A - Atf) ' { ) 

Fig. ^(a) shows the time evolution of u(t) obtained from a typical Lorenz attractor. Com- 
pared with the Rossler system (Fig. |l](b)), it clearly shows a broad distribution of u(t). 
Note that the scale of ordinates is different for Lorenz and Roossler system in Figs. |I|(a),(b) 
considering their difference of mean frequencies. 

For coupled Lorenz systems, steady states, synchronized oscillations |TD[, and front mo- 



tions [[LI]] have been studied. In particular, a linear stability analysis shows the existence 
of a large number of steady states |[L0|| . However it has not been clearly explained what 
characteristics of the system cause such dynamics. The purpose of the present work is to 
suggest the underlying physical mechanism and confirm its universality. 

While the synchronization of a small number of coupled Lorenz oscillators is achieved 
with weak coupling, this is not the case for a large number of coupled oscillators ||. In 
particular, a linear stability analysis shows that the Lyapunov exponent of the kth Fourier 
component is given by = A — 4e sin 2 (7r/c/iV) ||, where A corresponds to the Lyapunov 
exponent of the uncoupled single oscillator. Therefore for a fixed coupling constant e, as the 
system size iV increases beyond some critical value, synchronous states become unstable. 

Shown in Figs. 0(a) and (c) are numerical results of Eq. (0). Two-phase domains 



corresponding to the left lobe and the right form. In other words, the trajectory of each 
oscillator is confined either to the left lobe or the right in contrast to the uncoupled single 
oscillator. The domains show two types of dynamics: frozen or pulsating, as in Figs. 0(a), (c). 
In a frozen state, oscillations are suppressed to stationary states. However pulsating domains 
show recurrent oscillations [|22|| . Note that frozen and pulsating domains are coexisting. 

Now we will focus on the mechanism of such dynamics. As explained above, the Lorenz 
system has broad time scales. In Eq. ([!]), — 2e£j,— 2«/j and — 2ezi provides some additional 
damping, however, the forcing from the neighboring oscillators has minor effects since it is 
nonresonant due to the broad time scales. Therefore damping is dominant and oscillations 
are suppressed. 

Although we deal with the identical oscillators, these two types of domains can be well 
understood using the terminology of coupled nonidentical oscillators. It has been reported 
that coupled nonidentical oscillators with various frequencies |23||r^| exhibit very rich dy- 



namics including steady states such as oscillation death, as well as unsteady states such 
as periodic, quasiperiodic or chaotic evolution. Although we consider the identical Lorenz 
oscillators, they have broad time scales, so effectively they act as oscillators with various 
natural frequencies P0 |. We can then view the frozen states of coupled Lorenz systems as 



the oscillation death, and the pulsating states as the unsteady states of coupled oscillators 
of various frequencies. 

In coupled periodic oscillators p3) or chaotic oscillators such as Rossler systems |L2| , for 



the forcing from the neighbor to be nonresonant so that the oscillation death occurs, coupling 
with nonidentical oscillator is required. However, for the Lorenz systems, oscillation death 
can be obtained from the coupling of the identical oscillators since they have broad time 
scales. 

For a more quantitative study, the local frequency gradient, 5w = {\uii — u)i-i\ + \uii — 
Ui + x\)/2, is calculated from the same data of Fig. 0(a), and plotted in Fig. 0(b). We want to 
emphasize that uo is defined in the absence of coupling to reflect the intrinsic property of the 
local dynamics and given by Eq. @) which is derived from the equations of the uncoupled 



single Lorenz oscillator. Compare Fig. 0(a) with Fig. 0(b), and see Fig. 0(d) and (e), 
where u>i are obtained for the regions denoted by "frozen" and "pulsating" in Fig. 0(a), 
respectively. From these results, one can see that the oscillation death corresponds with 
large frequency mismatches between neighboring sites, and oscillating state with the small 
mismatches. Note that small clusters with close frequencies in Fig. 2(d) correspond to the 
boundary of domains where simple comparisons of frequencies between neighboring sites are 
not sufficient for the analysis of dynamics. 

The relation between the dynamics of domains and the frequency mismatches is more 
clearly shown in Fig. 0, where probability distribution function of frequency mismatches, 
Sw, of space-time points for frozen and pulsating domains are denoted as solid and dashed 
lines respectively To get the histograms in Fig. [| 100 independent numerical simulations 
with different initial conditions are considered. Initial transients (t < 300) are excluded, and 
data are obtained with a time step At = 0.025 for each space point during next 50 sec. Also, 
boundary points of domains are excluded. Fig. |3] shows that the probability distribution 
function for the pulsating domains is peaked at lower values of frequency mismatch compared 
with that of the frozen domains. 

The results are consistent with those of previous studies on the coupled nonidentical 



chaotic oscillators In that case, it was reported that the large frequency mismatches 
induce the oscillation death and the small ones bring the synchronized oscillation. 

Rigorously speaking, the exactly stationary states cannot coexist with the oscillating 
states in the same lattice. But the effect from neighboring pulsating domains are so well 
blocked that the frozen domains can be essentially stationary except at the boundary between 
the frozen and pulsating domains. For the frozen domains, we could not find any noticeable 
oscillations numerically. 

Also note that bistability of the local dynamics is important for oscillation death. Nu- 
merical results show that larger domains tend to be pulsating rather than frozen. It means 
that there is an upper limit for the width of the frozen domain. Without bistability, there 
could be no domain structure, and so oscillation death could not occur. 



Next, we investigate the effects of forcing on the coupled Lorenz systems, which confirm 
the above explained role of broad time scales again. For the preservation of the mirror 
symmetry of the system, a force Acos(Qt) is applied to the z variable. For uncoupled 
systems, synchronization with the external periodic forcing is imperfect due to the broad 
distribution of the intrinsic time scale [2D||. Here, synchronization is defined by the absent 



or negligible difference between the external forcing frequency Q and the average phase 
velocity uJ of the system obtained from a sufficiently long time series. For a fixed A, UJ — Q 
is calculated as a function of Q; typical results are presented in Fig. §(a). There exists a 
horizontal plateau with some fluctuations near zero. This imperfect synchronization region 
as a function of the external forcing frequency Q and amplitude A is depicted by solid lines 
in Fig. |(b). 

For coupled systems with time- varying forcing, clearly, oscillation death is impossible. So 
we concentrate on the effects of forcing on the confinement. The external force entrains the 
system, making a characteristic time scale. Hence the chaos suppression effect, which was 
enhanced by the broad time scales, is weakened. Thus the resonant external force breaks 
the confinement. As a function of uo and A, two regions (I and II) are obtained in Fig. ^(b). 
In Region I (shaded region), confinement is broken and the dynamics of the system becomes 
turbulent (Figs. U(d),(f)); In Region II (unshaded region), one can see two-phase domains 
with confined dynamics (Fig. f|(e)). Region I coincides with the region of synchronization 
with external forcing of the uncoupled oscillator. 

However there are some deviations: Region I is shifted to higher frequencies compared 
with the synchronization of the local dynamics. In fact, the situation is very similar to the 
results from the complex Ginzburg Landau equation with external forcing where coupling 
changes the frequency of the system, so that Arnold tongue structure changes compared with 
the uncoupled cases P^| . The UPOs of the Lorenz system typically consist of fast rotations 
in one lobe and slow transition to another lobe |2(J . With coupling, the long time transition 



between two lobes is prohibited and effectively the frequency of system is increased, which 
may explain the above described deviation. When the confinement is broken by the external 
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force, the prevented transitions are possible and the dynamics of the system becomes slow. 
Then the system can be nonresonant to the external force and confined again. Thus, at the 
boundary between Region I and Region II, one can see mixed patches of domain structure 
where domains with the confinement and without the confinement are coexisting as in Fig. 
1(f)- 

Thus far, we have shown that the broad time scales induce oscillation death and confine- 
ment in the coupled Lorenz oscillators. To confirm the universality of these results, we now 



turn to the dynamics of coupled Chua [gfj and periodically forced Duffing oscillators. 
For the coupled Chua oscillators, we consider the dynamics of the following equations: 

^ = c x {yi -Xi- g(xi)) + e(x i+ i + x { -\ - 2xi) (3) 

^ = c 2 {xi -yi + Zi) + e(y i+ i + yi-i - 2yi) 

^ = -CaVi + e(z*+i + - 2z { ), 

where g(x) = m\X + m °~ mi (\x + 1| — \x — 1|) (ci = 15.6,c 2 = l,m = —8/7, mi = — 5/7, c 3 = 
24,e = 5.0). 

The single Chua oscillator without coupling yields a chaotic attractor with two wing 
structure as in Fig. |5|(a). To check the broadness of the time scales, we should define the 
phase and calculate the instantaneous frequencies. However, in this case, a unique center of 
rotation is not well defined, so a simple formula for the instantaneous frequency as a function 
of arbitrary x, y, z can not be obtained. We define the phase for the Chua oscillator by the 
empirical- mode-decomposition method p5| , and then the instantaneous phase velocity shows 
broad time scales as in Fig. |5|(b). The dynamics of the coupled Chua oscillators is quite 
similar to that of the coupled Lorenz oscillators. For weak coupling, we have imperfect 
confinement and no oscillation death as in Fig. ^|(a). However, with increasing forcing 
strengths, confined domains show two types of motions : oscillation death or nearly periodic 
oscillations as in the case of Lorenz oscillators. As seen in Figs. |6|(b) and (c), frozen domains 
and oscillating domains coexist. 

The equations for the lattice of forced Duffing oscillators are given by 



dxi 

~dt 

dyi 
dt 



Hi + e(x i+1 + Xi-x 



ax\ + bxi + cyi + f sin(ut) + e(y i+1 + y^i - 2yi) 



2xi) 



(4) 



. The orbit of a single forced Duffing oscillator shows a chaotic bistable nature as in Figs. 
0(a), (c). But the situation is different with the Lorenz oscillator or the Chua oscillator. 
For the lattice of forced Duffing oscillators, oscillation death cannot be obtained due to 
the external time-varying force. Also, confinement is expected to be easily broken since 
the time scales of dynamics are dominated by the external forcing frequency. Since we are 
interested only in the confinement phenomena, we use the binary representation of space time 
evolutions of dynamics: white(black) denotes positive (negative) x in Figs. (b),(d),(e),(f). 
For weak forcing amplitudes, confinement is obtained with weak couplings (See Fig. 0(b)). 
But boundaries between the domains are oscillating. With increasing coupling constant e, 
width of domains increase, however the boundary instability does not disappear. For strong 
forcing amplitudes, it is difficult to obtain the confinement. For lower coupling strengths, 
spatiotemporal chaotic pattern occurs as in Fig. 0(d). With increasing diffusive coupling 
constant e (about more than e = 50), we get confined domains, however, confinement is 
still imperfect at the boundaries of domains as in Fig. 0(e). Further increasing coupling 
strength, the domain structure disappears and the confinement is broken as in Fig. 0(f). 
For sufficiently large coupling strengths, synchronization of whole system occurs. 

In summary, we have investigated the mechanism of confinement and oscillation death 
in coupled chaotic bistable systems. With the broad time scales of the local dynamics, 
diffusive coupling suppresses the oscillation and causes confinement and oscillation death. 
The result implies that in coupled chaotic oscillators, the general properties of dynamics 
can be classified by the characteristics of the individual oscillators. The author wishes to 
thank H.T.Moon, S.-O.Jeong and T.-W.Ko for many stimulating discussions and M. Baer 
for useful comments. This work was supported in part by the interdisciplinary research 
(Grant No. 1999-2-1 f 2-002-5) program of the KOSEF and in part by the BK21 program of 
the Ministry of Education in Korea. 



S 



REFERENCES 



[1] A.T.Winfree, J.Theor.Biol. 16, 15 (1967). 

[2] Y. Kuramoto, Chemical Oscillations, Waves and Turbulence (Springer- Verlag, 1984). 

[3] A. S. Mikhailov, Foundations of Synergetics I (Springer- Verlag, 1990). 

[4] K. Kaneko, Phys. Rev. Lett. 63,219 (1990); ibid. 65,1391 (1990); ibid. 69,905 (1992). 

[5] Gang. H and Zhilin. Q, Phys. Rev. Lett. 72,68 (1994); R. O. Grigoriev, M. C. Cross, 
and H. G. Schuster, Phys. Rev. Lett. 79,2795 (1997); N. Parekh, S. Parthasarathy, and 
S. Sinha, Phys. Rev. Lett. 81, 1401 (1998); N. F. Rulkov, Phys. Rev. Lett. 86, 183 
(2001). 

[6] A. Lemaitre, H. Chate, and P. Manneville, Phys. Rev. Lett. 77, 486 (1996); A. Lemaitre 
and H. Chate, Phys. Rev. Lett. 82, 1140 (1999). 

[7] L. Brunnet, H. Chate, and P. Manneville, Physica (Amsterdam). 78D, 141 (1994); R. 
Hureta, M. Bazhenov, and I. Rabinovich, Europhys. Lett. 43, 719 (1998). 

[8] J. F. Heagy, L. M. Pecora, and T. L. Carroll, Phys. Rev. Lett. 74, 4185 (1995); Phys. 
Rev. E 50, 1874 (1994). 

[9] A. Goryachev and R. Kapral, Phys. Rev. Lett. 76, 1619 (1996); A. Goryachev, H. Chate, 
and R. Kapral, Phys. Rev. Lett 80, 873 (1998); ibid. 83, 1878 (1999). 

[10] N. J. Balmforth, C. Pasquero, and A. Provenzale, Physica (Amsterdam). 138D, 1 
(2000); K. Josic and C. E. Wayne, J. Stat. Phys. 98, 1 (2000). 

[11] D. Pazo, N. Montejo, and V. Perez-Munuzuri, Phys. Rev. E 63, 066206 (2001). 

[12] G. V. Osipov et al, Phys. Rev. E. 55, 2353 (1997); ibid 58, 7198 (1998). 

[13] L.M.Pecora and T.L.Carroll, Phys. Rev. Lett. 64, 821 (1990); Phys.Rev.A 44, 2374 
(1991). 



9 



[14] M. G. Rosenblum, A. S. Pikovsky, and J. Kurths, Phys. Rev. Lett. 76, 1804 (1996); 
ibid. 78, 4193 (1997) ; K. J. Lee, Y. Kwak, and T. K. Lim, Phys. Rev. Lett. 81, 321 
(1998) ; E. Rosa, E. Ott, and M. H. Hess, Phys. Rev. Lett. 80, 1642 (1998). 

[15] V. S. Anishchenko, A. N. Silchenko, and I. A. Khovanov, Phys. Rev. E. 57, 316 (1998). 

[16] J. S. Park and K. J. Lee, Phys. Rev. Lett. 83, 5393 (1999). 

[17] S. H. Strogatz Nonlinear dynamics and chaos (Addison Wesley, 1994); E. Ott Chaos in 
dynamical systems (Cambridge, 1993). 

[18] A. S. Pikovsky et al, Physica (Amsterdam). 104D, 219 (1997). 

[19] E. N. Lorenz, J. Atmos. Sci. 20, 130 (1963). 

[20] M. A. Zaks et al, Phys. Rev. Lett. 82, 4228 (1999); E. H. Park, M. A. Zaks, and J. 
Kurths, Phys. Rev. E. 60, 6627 (1999). 

[21] D. Auerbach et al, Phys. Rev. Lett. 58, 2387 (1987); C. Grebogi, E. Ott, and J. A. 
Yorke, Phys. Rev. A 36, 3522 (1987). 

[22] H. T. Moon, Phys. Rev. Lett. 64, 412 (1990). 

[23] R. E. Mirollo and S. H. Strogatz, J. Stat. Phys. 60, 245 (1990); P. C. Matthews and 
S. H. Strogatz, Phys. Rev. Lett. 65, 1701 (1990); J.T. Ariaratnam and S. H. Strogatz, 
Phys. Rev. Lett 86, 4278 (2001). 

[24] H.-K. Park, Phys. Rev. Lett. 86, 1130 (2001). 

[25] T. Yalcinkaya and Y.-C. Lai, Phys. Rev. Lett. 79, 3885 (1997); N. Huang el al, Proc. 
R. Soc. London A 454, 903 (1998). 

[26] T. Matsumoto, L. Chua, and M. Komuro, IEEE Trans. Circuits Syst. 32, 798 (1985). 



10 



FIGURES 

FIG. 1. Instantaneous frequencies uo{t) obtained from (a) Lorenz system and (b) Rossler 
system. For the Rossler system, the equation is given by x = —y — z,y = x + ay, z = xz — cz + b 
(a = 0.2, b = 0.4, and c = 8.5) and the phase is defined by tan _1 (y/x). 

FIG. 2. Coupled Lorenz system, e = 1.0,7V = 300. Time evolution of x,- L is plotted in (a) 
(from i = 1 to i = 300) and in (c) (from i = 201 to i = 300). In (b), local frequency gradient 
5w = — + \uJi — Wj+i|)/2 is calculated from the same data in (a). In (a) and (b), the 
brightness is proportional to the largeness of u and Sw. In (d) and (e), u>i is obtained from the 
regions denoted by "frozen" and "pulsating" in (a): (d) frozen (e) pulsating. 

FIG. 3. Probability distribution function of the local frequency 

gradient, Sw = (\u>i — + \a>i — Wj + i|)/2, for frozen (solid line) and pulsating (dashed line) 

domains. 

FIG. 4. Forced Lorenz system, e = 1.0 (a) ZJ — Q vs Q, (b) Inside the solid line, imperfect 
synchronization occurs for the uncoupled single oscillator. In Region I (shaded region), the con- 
finement is broken and dynamics like (d) or (f) occur. In Region II (unshaded region), confined 
two-phase domain structure such as (e) forms, (c) A closer view of (b): triangle - (d), square - (f), 
elsewhere - (e). In (d),(e) and (f), the brightness is proportional to X{. 

FIG. 5. (a) Chaotic orbit of a single Chua oscillator and (b) its instantaneous frequencies u(t). 

FIG. 6. Coupled Chua oscillators (N = 300). Space-time diagram : (a) e = 1.0 and (b) 
e = 5.0, the brightness is proportional to x«. (c) Time series of Xj at i = 20, and % = 90 in (b) 

FIG. 7. Lattices of forced Duffing oscillators (a = -l,b = l,c = -0.25, u = 1, andiV = 300). 
Orbits of uncoupled Duffing oscillators for (a) / = 0.29 and (c) / = 0.40. Binary representations 
of space-time evolution for the coupled lattice are given for / = 0.29 ((b) e = 1) and for / = 0.40 
((d) e = 5, (e) e = 60, and (f) e = 100). 
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